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Testing the chiral behavior of the hadron spectrum* 
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'^Physics Department, University of Washington, Seattle, WA 98195 U. S. A. 

We analyze the chiral behavior of the hadron spectrum obtained with quenched Wilson fermions on 170 32'^ x 64 
lattices at /3 = 6.0. We calculate masses of hadrons composed of both degenerate and non-degenerate quarks. 
We reduce the statistical errors in mass splittings by directly fitting to the ratio of correlation functions. We 
find significant deviations from a linear dependence on the quark mass, deviations that are consistent with the 
higher order terms predicted by quenched chiral perturbation theory. Including these corrections yields splittings 
in the baryon octet that agree with those observed experimentally. Smaller higher order terms are also present in 
rUp and mjv. By contrast, the decuplet baryons are well described by a linear mass term. We find the decuplet 
splittings to be 30% smaller than experiment. We extrapolate our data to a by combining with the GFll 
results, and the best fit suggests that the quenched approximation is only good to 10 — 15%. 



1. TECHNICAL DETAILS 

We calculate quark propagators using the sim- 
ple Wilson action and both smeared Wuppertal 
and Wall sources. From these we construct three 
types of hadron correlators: wall source and point 
sink (WL), Wuppertal source and point sink (SL), 
and Wuppertal source and sink (SS). We use five 
values of quark mass given by k = 0.135 (C), 
0.153 {S), 0.155 (C/i), 0.1558 (C/2), and 0.1563 
(1/3), corresponding to pseudoscalar mesons of 
mass 2835, 983, 690, 545 and 431 MeV respec- 
tively where we have used 1/a = 2.33 GeV for 
the lattice scale. We use the three light quarks to 
extrapolate the data to the physical isospin sym- 
metric light quark mass m = (m„ -I- TO(j)/2. The 
physical value of strange quark lies between S and 
Ui and we use these two points to interpolate to 
it, and use C for the charm mass. 

Our overriding fitting criterion is to include 
as many time-slices as possible. We have not 
succeeded in developing an automated procedure 
that meets this objective and works in all cases 
when using two mass fits or incorporating the 
full covariance matrix. For the sake of unifor- 



Figure 1. Gomparison of mcs{t) for UiUiUi nu- 
cleon correlators with SL and WL sources. 
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* Poster presented by Tanmoy Bhattacharya and Rajan 
Gupta. These calculations have been done on the CMS 
at LANL as part of the DOE HPCC Grand Challenge 
program, and at NCSA under a Metacenter allocation. 



mity, therefore, we use results from single mass 
fits keeping only diagonal elements of the corre- 
lation matrix. In the cases we have checked, this 
does not significantly effect either the results or 
error estimates. All errors are obtained using sin- 
gle elimination Jackknife. 

We find that masses extracted from WL cor- 
relators are systematically lower than those from 
SL or SS correlators. For the pion and rho chan- 
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nels, this difference is ~ la, while for the nucle- 
ons it is ^ 2 — 3fT as exemphfied in Fig. This 
difference arises because the signal dies out be- 
fore contamination from excited state is fully re- 
moved. In cases where the signal persists to large 
enough times, the estimates from the wall and 
Wuppertal sources converge to a value roughly in 
the middle. Since SL and SS results are highly 
correlated, we take as our best estimate the aver- 
age given by {2xWL + SL + SS)/4. 

To extrapolate the hadron masses towards the 
chiral limit, and to test the forms predicted by 
quenched chiral perturbation theory (CPT), we 
make fits to both 1/2k and the non-perturbative 
definition of lattice quark mass described in 
Ref. The two give virtually identical results. 

2. MESON SPECTRUM 

A linear fit to using the six lightest {UiUj} 
points gives M2 = 0.0013(5)-h2.296(ll)mg. With 
our current data we cannot determine whether 
the non-zero (though tiny) intercept is due to fi- 
nite size effects, quenched chiral logs or the ab- 
sence of chiral symmetry with Wilson fermions. 
The linear fit to Mp (0.3296(59) -I- 2.54(14)™,) is 
shown in ^. The data show clear curvature, so we 
fit adding two types of corrections: mP^"^ (non- 
analytic term due to chiral loops) and m^. Both 
succeed in fitting the {UiUj} and {SS, SUj} 
points. These fits give 

Kc = 0.157131(9), 

= 0.157046(9), 
a-i(Mp) = 2.330(41) GeV (m fit) , 
a-i(Mp) = 2.365(48) GeV {m^/^ fit) , 
a-i(Mp) = 2.344(42) GeV (m^ fit) , 

Since the three fits to Alp have comparable x^, 
we cannot study CPT in detail. 

Using linear fits for pseudoscalar and vector 
states we determine the lattice strange quark 
mass by first extrapolating Mj^/M^, Mk* /Mp, or 
M^/Mp to m, and then linearly interpolating be- 
tween Ui and S until these quantities match their 
physical values. The results for {MS scheme, 
run down to 2 GeV) are ms{MK) = 89(2) MeV 



Figure 2. A linear fit to lightest six {UiUj} points. 
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versus m,{M^) = 106(6) MeV, i.e. a 20% dif- 
ference. Using M'^/M^ to fix uis implies that 
TOs = 25rn as we use the lowest order chiral ex- 
pansion to fit the data. On the other hand, using 
M^/Mp gives ms/m « 30, in surprisingly good 
agreement with the next-to-leading chiral result 
The differences presumably result from a 
combination of quenching and discretization er- 
rors. We use ■ms{M^) here. 

3. BARYON SPECTRUM 

We calculate three types of correlation func- 
tions made up of fiavors A,B,C. The spin-1/2 
S type (which also includes S and N) are la- 
beled A{BC} (symmetric in B, C); the A-like are 
A[i?C] (anti-symmetric in B,C); and the spin- 
3/2 decuplet are {ABC}. With 4 flavors S,U^, 
one can write 40 correlators of types A{BC{ and 
A[BC\ each. For B ^ C, SU(2) symmetry is 
broken and there is mixing between the A and E 
states. However, as explained in Q], for 5mt << 1 
(as is true for our data) the mixing can be ignored. 

To study mass splittings we make fits to ap- 
propriate ratios of correlators for a given source 
and sink. For example, fits to VY.{t) /T Nit) ~ 
^-{M^-MM)t^ give (Ms - Mn) directly. This 
has the advantage of both reducing some of the 
sources of systematic errors and of improving the 
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statistical errors (the errors in AM are 3 — 5 times 
smaller than those obtained from individual fits) . 

Our fits are motivated by the results of 
quenched CPT worked out in Ref. Q and are 
labeled chiral (non-analytic rr?J^ terms) and an- 
alytic (m^ terms). We first consider the octet 
hyperfine splitting, "E — A" , given by 
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BB 



(-8i:>/3) +ci 



+ di{mA + mB) + d[mB , (1) 



where Mab is the mass of the meson with fla- 
vor AB, etc. The constants Ci and di can be 
expressed in terms of parameters of the quenched 
chiral Lagrangian. For reasonable choices of these 
parameters, one expects \c[\ ^ |ci|. There is no 
useful information concerning di or d'^. 

Wc fit our results in two ways. First, we as- 
sume that the ci term is dominant, and make 
a linear fit (called chiral) to the 12 SU(2) sym- 
metric non-degenerate points as shown in Fig. ^ 
The sizable non-zero slope shows that significant 
terms of higher order than linear in the quark 
mass are needed. In our second fit we test to see 
whether our data can be represented as well with 
analytic terms. We find, by trial and error, that 
Eq. |l| with ci = c'l = and di = d'l works, as 
shown in Fig. |[ The data show significant cur- 
vature, so we have extrapolated to the physical 



point using a quadratic fit, corresponding to 
terms in baryon masses! The extrapolated value 
of Ms — Ma from the two fits are given in Ta- 
ble 1^; these agree with each other and with the 
experimental splitting. 

Next we consider the "S — iV" splitting. 
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4 ""'r r + d2{mA + ms) + d'^mc . (2) 



CPT suggests that 14 



< C2 



Thus we fit 



the data assuming C2 is the dominant coefhcient 
(chiral). Our best trial and error attempt with 
analytic corrections assumes C2 = c'2 = and 
d2 — ^2/4. Again there is definite evidence for 
curvature even though the collapse of data on to 
a single curve is not as good as in the S — A case. 
The difference in estimates of Ms — M^r given in 
Table |l| is indicative of this. 

Thirdly, we consider the difference "S — N" : 
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Here there is no expectation that C3 and C3 should 
be substantially different in magnitude. Never- 
theless, our "chiral" fit assumes C3 is the domi- 
nant coefficient, and fits the data reasonably well. 
A quadratic fit to the average quark mass (i.e. 
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Table 1. Octet mass splittings. Our most reliable 



Fit 


"m3/2" 


"m^ " 


Expt. 


Ms - Mn 


223(15) 


269(22) 


253 


Me - Mn 


339(20) 


358(21) 


375 


Ms - Ma 


78(14) 


84(9) 


77 



Figure 5. Linear fit to spin-3/2 baryon data. The 
four degenerate cases are shown with octagons. 
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assuming = d'^) is slightly better. Again, 
the extrapolated values in Table |l| are consistent. 
The important point is that in all three cases the 
higher order corrections are substantial. 

Since keeping only terms linear in the quark 
masses provides a poor description of most mass 
differences, we include an rriq^^ (or m^) term in 
our fits to the four degenerate combinations yield- 
ing Mn = 1070(35) MeV (1072(31) MeV). In 
each case we have included the same form for the 
higher order corrections in Mp when determining 
a and to. For our best estimate we take the mean. 
Mat = 1071(35) MeV. 

In the spin-3/2 baryons we can form 20 states 
with our four masses, and none of these mix with 
each other. As shown in Fig. ||, the linear fit 
Mao = 0.581(13) + 1.13(5)(toa + ms + mc)a 
works for all 20 quark combinations ABC. The 
mass differences calculated from this fit turn out 
to be w 25% smaller than the experimental val- 



ues. Note that this discrepancy is worsened if we 
use TOs(Mif). Linear extrapolation of Ma — Mn 
calculated from the ratios of correlators gives 
318(30) MeV. Including an m,^^ or to^ term in 
both Ma - Mn and Mp fits gives 365(44) and 
347(39) respectively. These estimates are slightly 
higher than the experimental value 293 MeV. 

4. CONTINUUM RESULTS 

Our best estimates at /9 = 6.0 with linear ex- 
trapolations in TOg using only the Ui quarks are 



Mn/Mp = 1.412(35) 
Ma/Mp = 1.800(47) 
Ma/Mn = 1.275(36) 



Expt: 1.22, 
Expt: 1.60, 
Expt: 1.31. 



The GFll collaboration has claimed, based 
on data aX (3 — 5.7, 5.93, and 6.17, that these 
ratios have a significant slope when extrapolated 
to a = 0. We update their "012" and "4" sink 
fits. The change, on adding our point, to the 
"012" fit is from Mn/Mp = 1.28(7) to 1.30(6) and 
from Ma/Mp = 1.61(8) to 1.62(7). The xVrfo/ 
for the new fits are 2.1 and 0.85 respectively. 
The analogous numbers for the sink "4" data are 
1.33(9) ^ 1.38(7) and 1.68(10) ^ 1.73(10) with 
X^/dof equal to 1.2 and 0.86 respectively. 

The difference between the two fits is due to the 
"012" and "4" sink data at /3 = 5.7. If we neglect 
the point at strongest coupling, (3 — 5.7, then the 
remaining three points again show smaller a de- 
pendence, similar to the sink "4" case. Thus, our 
preferred estimates are Mn/Mp — 1.38(7) and 
Ma/Mp = 1.73(10) from sink "4" fit. The am- 
biguity in the extrapolation makes it clear that 
data at more values of (3 are needed in order to 
reliably determine continuum results. 
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